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Abstract 

We study thermodynamic properties, specific heat and susceptibility, of 
XY quantum chains with coupling constants following arbitrary substitution 
sequences. Generalizing an exact renormalization group transformation, orig- 
inally formulated for Ising quantum chains, we obtain exact relevance criteria 
of Harris-Luck type for this class of models. For two-letter substitution rules, 
a detailed classification is given of sequences leading to irrelevant, marginal or 
relevant aperiodic modulations. We find that the relevance of the same aperi- 
odic sequence of couplings in general will be different for XY and Ising quan- 
tum chains. By our method, continuously varying critical exponents may be 
calculated exactly for arbitrary (two-letter) substitution rules with marginal 
aperiodicity. A number of examples are given, including the period-doubling, 
three-folding and precious mean chains. We also discuss extensions of the 
renormalization approach to a special class of long-range correlated random 
chains, generated by random substitutions. 

1 Introduction 

Phase transitions and critical phenomena in Ising spin systems with (dis-)order of 
various nature (random, quasiperiodic, self-similar etc.) have been an active re- 
search area for many years. The main questions deal with the relevance of these 
kinds of disorder to the thermodynamics of different models and the characteriza- 
tion of new, disorder-induced universality classes. For random systems, the Harris 
criterion || gives an heuristic scaling argument for the relevance of disorder. This 
has later been generalized by Luck |24| to general aperiodic disorder. The argument 
is perturbative in nature and is expected to hold for weak disorder: Comparing the 
local shift of the critical point due to aperiodic modulations to the distance from 
criticality, disorder should be relevant if the so-called wandering exponent u> (mea- 
suring the geometrical fluctuations of the mean coupling constant) exceeds a critical 
value @ 

w c = 1 - l/(Du) . (1) 

Here, v is the correlation length exponent of the unperturbed system and the model 
is disordered in D co-ordinate directions. 



1 



Especially since the discovery of quasicrystals in 1984, the effect of (determin- 
istic) aperiodicity on the thermodynamical properties of different models has been 
the topic of numerous, mostly numerical studies. Universal behaviour was found for 
most quasicrystalline systems, like Ising models on Penrose and Amman-Beenker 
lattices ||, |3^, ||(| and also in three dimensions pifl . Marginal scaling, on the 
other hand, has been observed for the surface roughness of two-dimensional qua- 
sicrystals [jlCj |j. Compare also the review for further references. Although al- 
most all results corroborate the Harris-Luck criterion, this is nevertheless somewhat 
more subtle for aperiodic or (more general) correlated disorder than for randomly 
disordered systems. This is because it is in general not the fluctuations of the 
coupling constants directly, but of some related microscopic parameters (reduced 
coupling constants) that should be considered. While the statistical fluctuations 
of uncorrelated random couplings will normally lead to a similar fluctuational be- 
haviour in the (a priori unknown) "reduced couplings" , this is no longer guarantied 
if the coupling constants are correlated or even distributed according to a deter- 
ministic rule. The solvable Ising models on the Labyrinth Q provide an example 
where the strong correlations among the coupling constants due to the "rapidity 
line parametrization" enforces Onsager universal behaviour even for (according to 
pi] ) relevant fluctuations of the mean coupling constant. 

Analytical results have so far been obtained only for a small number of systems. 
With the exception of the Labyrinth models, which are, when solvable, somewhat 
non-generic in their aperiodicity, these are one-dimensional free-fermion models, 
like tight-binding models or quantum chains. Most results rely moreover on a spe- 
cial choice of the aperiodic orderings (like the Fibonacci model) which makes them 
applicable to efficient trace-map methods or renormalization techniques derived 
therefrom |l6| |l| . Independently of trace map properties, the surface magnetization 
of aperiodic Ising quantum chains with constant transverse field has been calcu- 
lated exactly for certain substitution sequences p9| . Only recently, a decimation 
procedure in real space has been introduced p0| (again for particular substitution 
rules) that later could be generalized to obtain analytically the scaling properties of 
the whole class of Ising quantum chains with coupling constants following arbitrary 
substitution rules JIlJ . This lead to an analytical confirmation of Luck's relevance 
criterion for these models. 

In this article, after an introduction of the model in section 2, we show in section 
3 how the renormalization approach, as formulated for the Ising quantum chains, 
can be extended to aperiodic XY spin chains. It turns out that, for a given se- 
quence of couplings, the influence of the induced disorder may be different in the 
two models. Nevertheless, fluctuations turn out to be the basic concept for the 
demarcation of relevant from irrelevant disorder. However, the fluctuations of the 
sequence of coupling constants and of the induced sequence of reduced coupling 
constants, that determine the critical behaviour here, behave - in contrast to the 
randomly disordered case - in general differently for aperiodic order. This way, the 
Harris-Luck relevance criterion may be adapted to XY spin chains or, equivalcntly, 
to tight-binding models with aperiodic hopping. We calculate the scaling exponent 
of the central spectral gaps at criticality and derive therefrom in section 4 (following 
p2| ) the critical scaling behaviour of the specific heat and the zero-field suscepti- 
bility. Connections to localization properties of tight binding models are briefly 
mentioned. In section 5, we show how known results from trace-map approaches 
can be rederived, clarifying their origin in this broader context. As examples, we 
also give some new scaling exponents for different aperiodic chains with marginal 
disorder. In section 6, an extension of the renormalization approach to random 
substitution rules is proposed; finally we conclude with a short discussion. 
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2 The model 



The system we are concerned with here is defined by the following quantum Hamil- 
tonian: 



N 



ii y x (>>.;•: • >;<.) ■ ( 2 ) 

3=1 

The coupling constants e*' v are site-dependent and the operators crj' y denote Pauli's 
matrices acting on the jth site. Boundary conditions may be chosen periodical 
(cjv+i = ci) or free (en — 0). 

For a general set of coupling constants, this model is equivalent with a free 



fermion field 21 pi), the fermionic excitation energies A g satisfying the linear dif- 



ference equations 

A q ^ = eU^li + e^l (3) 
A,^> = e^a+sf^i. (4) 

Defining 

r$=$ ^-1=4?-! (5) 

if = 4f if-, - 4U , (6) 

these equations decouple into the eigenvalue problems of two independent tight bind- 
ing models with aperiodic hopping. This decoupling can also be carried through on 
the level of the spin chain Hamiltonian itself and has been used there to analyse 
XY-chains with random bonds ||. Difference operators of the kind ^) underly 
various physical models and may also be interpreted as a phononic model with vary- 
ing spring constants or the transition matrix of a one-dimensional random walk in 
an aperiodic environment. The Ising quantum chain with transverse magnetic field 
in its fermionic form also gives rise to a similar set of equations, the field variables 
replacing the e y couplings. In a renormalization scheme has been defined for 
the case of a uniform magnetic field (or, more generally, field variables depending 
on the neighbouring coupling constants), thereby effectively decoupling the degrees 
of freedom that finally enter the renormalization scheme. In our situation, however, 
the e y couplings will not be determined through their neighbourhood, but, together 
with the e x couplings, follow the aperiodic sequence that defines the model. 

For uniform or randomly distributed coupling constants, the XY- chain exhibits 
a zero temperature phase transition from an X- to an Y-fcrromagnetically ordered 
phase at [ln£ x ] ,j = [lne y ] at; . For the quasiperiodic Fibonacci sequence, a non- 
universal scaling law has been found for isotropic couplings (ef = e y ) with scaling 
exponents depending on the coupling constants . 

Here, the site-dependent coupling constants e*' v shall be chosen according to a 
sequence generated by an arbitrary substitution rule 

g : at — > Wi (7) 

with words Wi taken from an n-letter alphabet A with letters ax, . . . ,a n . In the 
following, wf shall denote the £th letter and # Q (wi) the total number of letters 
a in Wi. Since also (non-overlapping) pairs of consecuting letters in a word will 
play an important role, we further define # Q( g(wi) to give the number of pairs (a/3) 
contained in Wi. Finally, let \wi\ be the length of W{. 

The parameters of the aperiodic model are given by the ratios of the £^ y 's. We 
define, on a logarithmic scale, 

ry=ln<+lne-ln<--lne (8) 
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and 



(9) 



The n-letter model thus contains n — 1 independent variables which parametrize the 
strength of isotropic aperiodicity and n parameters Aj which determine the aperiodic 
anisotropy of the model. For notational clarity, we will restrict discussions from now 
on to substitution systems on a two-letter alphabet, and comment only briefly on 
extensions to general ro-letter substitutions, which in most cases can be dealt with 
along the same lines. We assume the following normal form for the substitution 
rule 



Q ■ 



"'a 
U'b 



= aw a 
= bw' h 



(10) 



which does not mean a restriction as any two-letter substitution rule may be trans- 
formed accordingly without changing the model |p"lf . In principle, the special form 
for g is not needed to make our renormalization group work, but it simplifies some 
of the calculations. The two-letter model contains three independent parameters 



In 



£ £ u 

Ft — /I 



'b°b 



and 



A a 
A h 



He x Jel) 

H £ b/4) 



(ii) 



(12) 



Some basic properties of the sequence generated by g are already contained in the 
corresponding substitution matrix 



#a(w a ) 

#b{w a ) 



#a(Wfo) 

#b{w b ) 



(13) 



with eigenvalues A±. The leading eigenvalue, A+, gives the asymptotic scaling fac- 
tor of the system size with the number of iterated substitutions, the entries of the 
corresponding (statistically normalized) eigenvector i>+ determine the frequencies 
Pa : b of the letters a, b in the limit word |28| . The further eigenvalues of the substitu- 
tion matrix are connected to fluctuation modes present in the sequence. Especially, 
the next to leading eigenvalue (here A_ , of course) determines the fluctuation (or 
wandering) exponent 



In I A. 
In A i 



(14) 



For two-letter substitution rules in particular, a classification of the substitution 
matrices according to the fixed points of their reductions modulo 2 will be helpful. 
Up to the exchange of the letters a and b, there are five possible fixed points of 
[7W e ] moc j2, grouping into three cases, for which we will show the following: 



1. 



a) [M e ] mod2 = 

b ) [M e } mod2 







1 
1 



1 1 
1 1 



1 




(15) 
(16) 
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The chain is in its critical phase at T = iff 



Pa(#ab ~ #6a)(w a ) +Pb(#ab - #ba)(wb) 



(17) 



This condition holds in particular if (A Q # a +Af,#h)(w Q ) = ±r(# Q (, — #& a )(w Q ) 
where we find unperturbed scaling behaviour [z — 1). Otherwise, the rele- 
vance of the aperiodicity to the critical scaling is identical to the one for Ising 
quantum chains with couplings following the same aperiodic sequence: We 
obtain z = 1, marginal scaling with non-universal (1 < z < oo), or random- 
like behaviour (z = oo) for |A_| lower, equal, or larger than 1 respectively. 
Finally, off criticality, there is a gap in the fermionic integrated density of 
states (IDOS) at the band center. 



2. 



a) [M g ] mod2 

b) [M B ] mod2 



1 1 



1 
1 



1 

1 



1 1 
1 



1 

1 1 



(18) 



(19) 



We obtain a critical ground state and a vanishing gap in the center of the 
IDOS iff 



v + ■ A = 



(20) 



In particular, this implies that the model is critical for any kind of isotropic 
aperiodicity (A a = A/, = 0). In this case, the relevance of the aperiodic 
modulation is determined by the renormalization group eigenvalue \ xx = 
I#a6 — #ba|(w a Wfc) and is completely independent of the relevance of the same 
aperiodic ordering in the couplings of the corresponding Ising quantum chain. 
The relevance of aperiodic anisotropy, on the other hand, again agrees with 
the Ising case. If isotropic and anisotropic aperiodicity are both marginal, 
we obtain a nonuniversal scaling exponent 1 < z < oo depending on two 
parameters. 



3. 



[A^£)]mod2 



1 





(21) 



The criticality condition is again given by (20). Isotropic aperiodicity is ir- 
relevant if #ab(w b ) = #ba(w b ) or if 2# bQ (w b ) = #b(w b ). Here, # b o(wb) gives 
the number of 6's in w b with w b = wibw2 and # a (wi) even. In all other cases, 
the induced disorder is relevant and z = oo. Again, the relevance of aperiodic 
anisotropy agrees with the Ising case. 

Please keep in mind that discussions may be concentrated on the five fixed points, 
since any substitution rule can be transformed into an appropriate one by taking a 
suitable power. This clearly does not change the limit chain. In the following, we 
show how to derive these relevance criteria within an exact renormalization scheme. 



3 The renormalization group 

In this section, we generalize the renormalization procedure introduced for the Ising 
quantum chain in [|llj . It relies on a decimation process found in |2(j] for particular 
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substitution rules and uses a special star-product technique, originally developed 
many years ago in the context of scattering theory (see p9[|). 

We introduce S transfer matrices for the decoupled sets of difference equations 

§§ 



V21 J W Wl-l 



fhk-1 ] = S 
V21 



k\l 



V2k 
V21-1 



(22) 



Since both systems are related under the exchange of e x and e y couplings, we 
may concentrate on the rj equations from now on. The S matrices transform by 
*-products like 



S, 



k\l 



Sk\k+i * s< 



fe+l|fe+2 



l\l 



l-k 

^fc £fc+i-l|fc+i 
i—1 



where k < I and the ^-product of two matrices is defined as 



(23) 



ci 



h 

di 



"2 
C2 



b 2 
d 2 



a 1 






d 2 



_L fbicia 2 bib 2 

d 1 a 2 \ cic 2 dib 2 C2 



The form of the elementary matrices follows from (|[ ||). For e| fe _ 1 = e 



(24) 



2 k 



and e\ M 



z y we obtain 



^fc|fc+l = SafJ-f — 



(25) 



In the corresponding problem of the aperiodic Ising quantum chain with constant 
transverse field h = 1, each S matrix contains just one coupling constant (e% in 



(25), setting all e y — > h = 1), enabeling an easy renormalization procedure re- 
versing the substitution steps by * multiplication of the corresponding S transfer 
matrices. With an S matrix depending on three consecuting coupling constants, 
the substitution rule has to be redefined appropriately in order to make an inverse 
transformation possible by taking * products. As the part of a coupling constant 
in the Ising model is just taken by the ratio of two consecuting coupling constants 
in the XY case, a substitution rule on pairs of letters (actually rather than triples) 
is needed here to make the renormalization procedure work. The Ising problem on 
n letters thus corresponds to a problem of dimension n 2 in the XY case - if such a 
substitution rule can be found at allQ. For the moment, let us concentrate on the 
cases 1 and 2 in the above classification, where a substitution rule of the desired 
form is easily constructed. Here, \w a \ + \w@\ is even for any a, f3 € {a, b} and we 
obtain a pair substitution 



02 ■ (a/3) -> w af} = w a wp 



(26) 



with substitution matrix 



M 2 = 



/ #aa(Waa) #aa(w a b) #aa{w b a) #aa(w bb )\ 

#ab(Waa) #a,b(w a b) #ab(w b a) #ab(w b b) 

#ba{w aa ) #ba{w a b) #ba{w b a) #ba(w b b) 

\ #bb{w aa ) #bb(w ab ) #bb{w ba ) #bb{Wbb) ) 



(27) 



1 The pair substitution needed here is entirely different from the one used to describe the Ising 
quantum chain with coupling constants depending on the two end-points (site-problem), considered 
in b3|. In the Ising site problem, the chain is divided into overlapping pairs and each coupling 
constant appears in two pairs (with each of its neighbours), but only in one pair here. This may 
lead to a much more pronounced change in the scaling behaviour, see below. 



G 



Let Xi be the eigenvalues of M2, Ai = Xpf being the Perron- Frobenius one. We 
now may adjoin each transfer matrix S a f3 7 to the pair (a(3) and define corresponding 
reduced coupling constants 



Ma/3 



ln^-ln^. 



(28) 



The renormalization transformation is now obtained by reversing the substitution 
procedure (p6[) 



Saf3 j = 



±(el/ey)exp p, a/3 



± exp jl al 3 



■x-o.fll/2 
i=l ° 3 ° 3 



2i+l 

a/3 



(29) 



(30) 



where 



1+1 



J a/3 



= 7. In (J29|) , we have introduced weights n a p as additional parame- 
ters in our renormalization group. These account for the fact that the renormaliza- 
tion blocks in general will be asymmetric (which causes n + and k~ to differ), may 
contain different coupling constants and also vary in length, which results in differ- 
ent local weights K a p and n a '/3'- Note that formally, in (p9j), only the e x couplings 
are renormalized, while the e^'s just keep their values. 

The renormalization flow on the critical surface (defined through the vanishing 
of the "mass gap" , A = 0) may now be directly obtained from (|24| ) as 



where, in the 2-letter case, 







f 2A a \ 


f-ab 


1 


A a + A b + r 


fJ-ba 


~ 2 


A a + A b - r 


\Vbb) 




V 2A b / 



(31) 



(32) 



The reduced coupling constants thus transform with the transpose of the pair substi- 
tution matrix Mi- Note that the RG transformations are the same in both decou- 
pled eigenvalue systems, however with different initial conditions (A a b — * — A a b ) of 
the reduced couplings. Obviously, the "Onsager fixed-point" of the uniform model 
just corresponds to fi = and the eigenvalues and eigenvectors of M.\ directly 
reveal the RG eigenvalues and scaling fields of the model. The contributions in 
the directions of different scaling fields are measured by the scalar products of \x 
with the eigenvectors ~V\ of Mi- In contrast to the renormalization of the Ising 
quantum chains the vector of reduced couplings fi is constrained here to a In — 1 
dimensional subspace of the vector space spanned by the frequencies of n 2 different 
pairs of n letters. It may thus happen that certain scaling fields vanish for arbitrary 
choice of fi which indicates that the corresponding fluctuational mode is not present 
in the problem. The contribution to the leading scaling field (with corresponding 
eigenvalue Xpp = A+, see below), however, does not vanish for a generical but 



/X • V PF = ^2p(a0)Hal3 = [ln£ X ] av ~ [ln£ y ] c 

(a/3) 







(33) 



leads to the well-known criticality condition for these models (e.g. [p7|). Any non- 
zero contribution immediately drives the system off the critical surface ((J, —> ±00) 
into the X or Y ferromagnetic phase. Note, however, that the whole model is critical 
if the criticality condition is fulfilled for just one of the decoupled subsystems. 
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The presence of aperiodic disorder in the chain leads to non-zero contributions 
in the direction of the additional scaling fields. Let A2 be the largest eigenvalue with 
a non-vanishing scaling field for critical couplings. If this scaling field is relevant, 
that is, iff the wandering exponent of the sequence of reduced couplings is positive, 

In I A 2 1 , . 

in Ai 

the system flows to the corresponding strong coupling fixed point of the RG, where 
the reduced couplings divide into two types, taking the values ±00 respectively. A 
(simple) eigenvalue | A2 1 = 1 leads to a marginal scaling field and the system flows 
to a fixed line with continuously varying exponents. This confirmes the Harris-Luck 
criterion for these models, since D = v = 1 for XY and Ising quantum chains, and 
thus uj c = according to ([!]). 

Before we discuss the effects of marginal or relevant aperiodicity to the critical 
behaviour, let us have a closer look at the spectra of the pair substitution matrices 
and the induced RG flows in the cases 1 and 2 above. It is worthwhile to consider 
for a moment the situation of vanishing isotropic aperiodicity, r — 0. Using 

2#qq(w Q)3 ) + #ab{w a p) + #ba{w a p) = # Q (W Q ) + #a{wp) (35) 

we recognize that the vector of the anisotropy parameters A transforms with the 
transpose of the original substitution matrix M B 

A = M l e A . (36) 

However, this means that the fixed point structure and the RG flow near the fixed 
points is identical for XY-chains with aperiodic anisotropy and the aperiodic Ising 
spin chain analysed previously [jlf). Of course, (|3rj ) implies that the spectrum of 
M g is contained in the one of Mi- These properties generalize without any change 
to n letter substitutions with \wi\ all even or all odd. Let us now see what happens 
to the RG flow when isotropic aperiodicity is turned on. 

• In case 1, where both \w a \ and are even at the fixed point of [Af e ] moc j2, 
the entries of Mi are related as 

# a p(w a >p>) = # a p(w a >) + # a p(wp>) ; a, (3, a', € {a, b} (37) 

and we conclude that the spectrum is just the set {A + ,A_,0, 0}. Using 
X+Pab — Pa#ab(wa) + Pb#ab(wb) we recognize that the criticality condition 
d33| ) reduces to equation (Jl^) in this case. Clearly, for substitution rules with 
# tt( 3(w 7 ) = #/3 Q (ui 7 ) Va,/3,7, isotropic aperiodicity is completely irrelevant. 
In any other case, isotropic aperiodicity and aperiodic anisotropy are coupled 
in ([IT]). The critical manifold (which always has to be symmetric under the 
exchange of indices x and y) splits into two submanifolds corresponding to the 
criticality of the decoupled subsystems which intersect at r = 0. Put in other 
words, isotropic aperiodicity deforms (and splits) the critical manifold, leading 
to a finite renormalization of the anisotropy parameters in the first RG step, 
but fixed point structure and RG flows remain unchanged otherwise. The 
changes of the critical behaviour induced by modulated aperiodicity remain 
closely related to the Ising case. This scenario generalizes without change to 
the n-letter case. 

• In case 2, |to | and \wb\ are both odd and we establish the following relations 

#ab(w aa ) = #ba{w a a)] #ab(w b b) = #ba(w b b) (38) 
#ab(w a b) - #ba(w a b) = #ba(Wba) - #ab(w b a) (39) 
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In absence of anisotropy this means that fi is already an eigenvector of M\ 
with eigenvalue X xx = (# a fc — #&a)(wa&), independent of the detailed form 
of the substitution rule. Note that X xx < A+, and \ xx = A + only in the 
degenerated case where w a b = (pb) n and the limit chain is periodic (with 
period ab). For all other cases, we conclude that the XX chain is critical 
for any amount of aperiodicity induced by substitution rules of the given 
form. In the XY case, the parameters r and A renormalize independently. 
In other words, the isotropic aperiodicity leaves the anisotropy parameters 
unrenormalized and does not deform the critical surface, but introduces an 
additional scaling field in the RG. Since this may be marginal or relevant, the 
scaling behaviour in general will be independent of the "Ising case". Note 
that the remaining eigenvalue A 4 = 5 \ Wn \ + 8 h„, , — 1, with eigenvector 

00 * w' a " ,a wl b ,b ' ° 

V4 = (1, —1, —1, 1)*, does not affect the RG transformation since \± ■ V4 = 0. 
Again, these properties generalize to n- letter substitutions with \wi\ all odd. 
Here, isotropic aperiodicity in general leads to contributions to n—1 additional 
scaling fields. 

So far, the third case in the above classification had been set aside. Here, \w a \ is odd 
and \wb\ is even at the fixed point of [M e ] mo d2- This makes things slightly more 
complicated since we cannot apply our pair substitution ( |2^ ) here. Nevertheless, an 
exact renormalization scheme may be set up also in this case, at least for two-letter 
substitution rules. The main idea is not to construct a substitution rule for pairs 
of letters, but for all substrings of the chain with an even number of a's and 6's 
and of minimal length (that is they cannot be divided into smaller strings with the 
same property). Obviously, (aa) and (bb) are examples for such minimal strings, a 
general string s with length 2k > 4 begins and ends with a pair of letters ab or ba 
with an arbitrary permutation of k — 1 pairs aa and bb in between: 




For a given substitution rule, the number of different minimal strings is always 
finite, hence a substitution rule on a finite "alphabet of different strings" Si can 
always be found 

g s : s t -> w Si = g(s})g(si) . . . (41) 

with substitution matrix 



[M s ] y = # s > Si ) (42) 

since w Bi may always be dissected into minimal strings. For real space renormaliza- 
tion, in a first step, we contract the strings by star multiplication of the correspond- 
ing S transfer matrices, this way assigning a single degree of freedom to each string. 
After that we proceed the usual way, reversing the substitution steps of the string 
substitution by decimation. Scaling fields and renormalization group eigenvalues are 
again determined by the action of the transpose of the substitution matrix M s on 
a scaling vector fi with entries corresponding to the different strings, where, from 
the initial conditions, 

fi St = A a # a (s t ) + A 6 # b (s 4 ) + r(# a6 - # 6a )(si) . (43) 

A more detailed analysis of the string substitution is given in the appendix, with the 
following results. As in the above cases, the anisotropy parameters transform with 



(|36[ ), leading to the Ising- like renormalization flow. For two- letter substitution rules, 
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the effect of isotropic aperiodicity is very similar to the case 2 above. A finite r does 
not lead to a renormalization of the anisotropy parameters, the criticality condition 
is given by (|2^). In particular, this means that the aperiodic XX model is always 
critical. But as in case 2, isotropic aperiodicity introduces a new scaling field in the 
model, with RG eigenvalue A s = 0, if (# a fc-#&a)(>&) = 0, and X s = (2# b0 ~# b )(w b ) 
otherwise. Again, the question arises whether these results generalize to the n-letter 
case. Note first that the method presented here can be applied only to a subclass of 
n-letter substitution rules. In general, a transformation into a string substitution 
with strings of even length will not be possible. Within these limits (but probably 
also beyond) the transformation of the anisotropy parameters A aj g generalizes just 
as in cases 1 and 2. On the other hand, isotropic aperiodicity in general will alter the 
criticality condition and renormalize A, but will always also introduce new scaling 
fields. 



3.1 Critical scaling behaviour 

The determination of the critical scaling behaviour of the lowest fermionic excita- 
tions may now be performed along the same lines as in the case of the aperiodic 
Ising quantum chain Jll]]. We will therefore only give a short account here. For 
cases 1 and 2, the RG transformations of the weights and fermion frequencies are 
to linear order in A 

A+ = M+A+ ; A_ = MA (44) 
where A± = A(k± q , k± , nf a , k±)* and 

kcsl/2 / y k-i \ 2 

k=l \ w 2 ^ 1 1 = 1 J 

\w a p\/2 /kval/2 \ 2 

M a0, a '/3' = £ ^f^.affi' II ^^w^w^) ' ( 46 ) 
k=l W=/s+l J 

For the derivation of (p5|), the special form of the substitution rule ( JToj ) is used. At 
the Onsager or marginal fixed points, a similarity transformation of M + yields the 
more symmetric form 

hajsl/2 k 

M %, a >f3> = exp(T2Ma/3) 5 ™l k f 1 K k l3 ><*'P' II ex P( ±2 ' i vgf x v%,') (47) 

k=i e=i 

In this form, the matrix elements are functions of the reduced couplings alone. The 
transformations in case 3 are analogous. The vectors of the weights and fermion 
frequencies converge under iteration of J44| ) to the Perron Frobenius eigenvectors of 
M and the scaling behaviour of the lowest fermionic excitations (q <C N) is then 
given by the normalization condition of A± as |]llf 



a f q Y ln(A M +A M -) 

A ^{n) 5 Z = 21nA + (48) 



where Xm± are the largest eigenvalues of Af ± . For irrelevant aperiodic modulations, 
we obtain 

M + = M = Ml (49) 
Am+ = ^M- = ^+ (50) 
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and hence z = 1. Since the fixed-point values of weights and couplings are indepen- 
dent of r and A and thus the same as for the uniform chain we may also conclude 
that the low energy excitations are equally spaced, in accordance with the predic- 
tions of conformal invariance. For the related problem of the Ising quantum chain, 
this has been observed numerically before (e.g. ||). 

Near the marginal fixed points, the coupling constants and thus the eigen- 
values Am± take non-trivial values. For aperiodic anisotropy (r = 0), the RG 
equations may exactly be reduced to the corresponding ones of the Ising quan- 
tum chain in a constant transverse magnetic field (h — 1), with coupling constants 
£ a b = ex P(Aa,fc)- Indeed, numerical observations show that not only the scaling 
exponents, but the entire low energy spectra are identical up to a common factor 
(altered fermion velocity). Note that this correspondence does of course not extend 
to the entire (high energy) spectrum. The scaling exponent z has been calculated 
explicitly for arbitrary two- letter substitutions in jllj] . After a proper renormaliza- 
tion of A, also the exponents of marginal substitutions of case 1 are found easily 
through this correspondence. For substitution rules of the second case, isotropic 
and anisotropic aperiodic modulation may be independently marginal. Generically, 
this leads to a scaling exponent 1 < z < oo depending on as many parameters as 
marginal scaling fields are present in the problem. Especially in the n-letter case, 
however, there are also exceptional substitutions, where marginal RG eigenvalues 
are in fact marginally irrelevant, leading to z — 1, or also lead to a marginally 
relevant scaling behaviour with z = oo, see jl2| for a discussion. 

For relevant aperiodic modulations, the reduced coupling constants do not tend 
to a finite limit, but (generically) grow with the second largest eigenvalue A 2 of M 2 - 
As a consequence, Xm± finally scales like Xm± ~ exp(cA2 ) in the nth renormaliza- 
tion step, resulting in a scaling behaviour of the lowest gaps as (see JTl] ] for a more 
detailed discussion in the Ising case) 

A q ~cM-c(N/qr*) (51) 



with u>2 defined in (34). Again, the same scaling behaviour, with the wandering 
exponent uj g of the sequence of couplings directly, had been found for the Ising 
quantum chains |2^| . Note that, contrary to random disorder, for aperiodic se- 
quences in general aj g ^ u>2 . 

In general, the RG flow to strong couplings may lead to rather unusual critical 
properties, where typical and mean values of various exponents (e.g. the correlation 
length critical exponent of the spin chain) no longer coincide H fL5|, see also the 
discussion below. The scaling exponent z to be calculated here is, however, not 
effected by "untypical events" of this type. 



4 Thermodynamical properties 

In this section the consequences of the scaling behaviour of the fermionic spectrum 
to the thermodynamics of the spin chain are discussed. This may be done in analogy 
to the analysis of the Fibonacci-XX-chain in (22] . 

The critical scaling of the low-energy spectrum directly implies the scaling form 
of the integrated density of states in the thermodynamical limit as |2^, [ll[] 

H(A) ~ A 1 / z 5 (lnA/lnA+) ; A^O (52) 

(g is a function with unit period) for marginal or irrelevant aperiodicity and 

H(A) ~ (ln|A|r 1/w2 ; A^O (53) 
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in the relevant case. The free energy (per spin) of the XY- chain at finite tempera- 
ture 1//3 is given by an integral transform of the fermionic IDOS as ^2| 

q 

= - J dH (A) ln(l + exp[/3A]) . (55) 



Now, the specific heat is given by 



B 2 B 2 r 

c * = ftp \-m = t J dH ^ 



A 2 



cosh'(/3A/2) 



(56) 



At low temperature, this expression is dominated by the small A region and the 
T — > scaling behaviour of C v is completely determined through the critical scaling 
of the fermionic spectrum 

C v ~T 1 l z G{\nT/\n\+) ; C v ~ l/(lnT) 1+1 /" 2 (57) 

for marginal (irrelevant) and relevant aperiodicity respectively, G is again a peri- 
odic function with unit period. Similarly, the susceptibility at vanishing field in 
z-direction may be derived to leading order as (with A(/i) = A(h = 0) + h • r(A), r 
bounded) 



d 2 f(h) 

Xz dh 2 



g f dH ^ ( 58) 
h=0 4 J cosh 2 (/3A/2) ' 



and 



^~T 1 / z - 1 G'(lnT/lnA+) resp. X * ~ l/(T[lnT] 1/w2 ) (59) 



Thus the susceptibility diverges for any marginal or relevant aperiodic perturbation. 
Note that for uj2 = 1/2, which is the mean fluctuation exponent for uncorrelated 
random disorder, these expressions coincide with the scaling behaviour of the ran- 
dom chain Q. 

As stated above, the aperiodic XY spin chain is essentially equivalent to two 
decoupled tight binding models with aperiodic hopping. In this context, the scaling 
exponent z calculated above determines the localization length at half filling. It 
is well known that the one-dimensional tight-binding model with random hopping 
exhibits a single delocalized state at the band center. This is also the case for any 
kind of aperiodic disorder (also due to random substitutions, see below) fulfilling 
the criticality condition. Using the Thouless relation [[33| we obtain from (^2|,^3|) a 
diverging localization length £\ at A — 0, like 

£ A ~ A" 1/z ; A (60) 

£\ ~ | In A|~ 1+1 / w ; A^O (61) 

for marginal (or irrelevant) and relevant perturbations respectively. A more de- 
tailed analysis of the resulting state (extended or critical) is possible for particular 
examples with the methods of jlTj. 



5 Examples 

• The Thue-Morse chain, generated by 

Qtm : { I Z fa ^ 
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belongs to case 1 and is critical whenever A a = —At,. The induced disorder 
is irrelevant, since # a b(w ab ) = #ba(w a b)- 



• The period doubling chain, generated by 



^ : ( b - 7a (63) 

is an example of a substitution rule of case 1 with |A_| = 1. The criticality 
condition reduces to 2A a + Ab + r = 0; note in particular that the pd XX chain 
is not critical. The critical scaling exponent is well known from the Ising case 
@ 

_ ln(2cosh(A a /2)) 
Zpd- (M) 

and is also for the XY- chain a function of only one variable. 

• The so-called precious mean chains jl3) are generated by substitution rules 
with a substitution matrix of the form 

M fc =(^'j) (65) 



n b 



with eigenvalues Afc± = (kzL^k 2 + 4)/2. According to the above classification, 
they all belong to case 2 and are critical for Afc_|_A a = — A/,. While criticality 
only depends on the anisotropy parameters, the critical exponent depends 
solely on r. Since |Afc_| < 1, anisotropic precious-mean modulations are 
irrelevant. On the other hand, it is straightforward to check that \\ xx \ — 
1, thus the isotropic aperiodicity is always marginal. For even k (fc = 2 
corresponds to the silver mean chain), we eliminate blocks corresponding to 
double substitution steps in the RG transformation and obtain the scaling 
exponent 



hi 9, „ 1 



Zk 



9 & = - + (66) 



where 



In A/. . 4 



p = exp(r/2) + exp(-r/2) . (67) 



For k odd, k = 21 — 1, in each RG step three substitution steps have to be 
reversed. We (finally) obtain the scaling exponent 



„ ln6i 1 
31nA fc+ 2 



(Pi(r)p 2 + /^w7+4) , (68) 



where 



= l 2 sinh[fr] + (l-l) 2 sinh[(l-l)r] 
sinh[r] 

and p as defined above. The first term of this series, with k = £ = Pi(r) = 
1, corresponding to the Fibonacci chain, had already been obtained in p2fl 
using the well-known properties of the Fibonacci trace map [^6|. Also for 
general precious mean chains, spectral scaling exponents may be calculated 
by_ trace maps due to the existence of invariants. This has been done in 
Note however that the scaling exponents found in Q (given in terms 



by t 
0- 
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of Chebyshev polynomials) do not simply translate to the above expressions 
since the transfer matrices of the aperiodic hopping problem do not have 
unit determinant. In contrast to the aperiodic potential problem [ fl4| , the 
scaling exponents here behave differently for k even or odd in the limit of weak 
incommensurability k — > oo. For fixed ratio of the couplings r, aperiodicity 
becomes irrelevant for k even (lim^oo Zk = 1), but not for k odd {zu — ► oo). 

The precious mean chains are just those quasicrystalline chains that result 
from the so-called cut-and-project formalism which the slope of the cut space 
given by Xk+ = [0, k]. By the sucessive application of different precious mean 
substitutions, a much more general class of cut-and-project chains may be 
generated. Indeed, since the eigenvector v xx to the marginal eigenvalue X xx 
of M2 is independent of k, the marginal scaling property holds also for this 
more general class of chains. Quadratic irrationalities in particular, which are 
observed in real quasicrystalline matter, are given by periodic continuous frac- 
tions and lead to cut-and-project chains that may be generated by a periodic 
application of precious mean substitutions. Thus also the scaling exponent z 
can be calculated using the above method. 

Finally, we like to stress that the origins of the marginal scaling behaviour 
observed for the interface roughness of Fibonacci surfaces |l(], ^ and the XY 
quantum chains on the other hand are independent. According to the Harris- 
Luck criterium, the former is connected to the fact that the unperturbed 
correlation length exponent which enters (|l|) is v — 1/2 there and the precious 
mean substitutions, being volume preserving, lead to the marginal wandering 
exponent u> e = — 1 fis}| . Substitution rules that lead to marginal scaling in 
only one of these situations are easily constructed; it is by chance that the 
precious mean chains fulfill both marginality conditions. 

Different types of the three-folding chain are defined by substitution rules with 
substitution matrix 

M3/ = (^) (70) 

This is one of the simplest examples of a chain where both isotropical and 
anisotropical aperiodicity are independently marginal (with exception of the 
special form g : a — > aba; b — > bab which leads to a periodic chain and z = 1). 
At criticality (A a = —A;, = A) we thus obtain continuously varying scaling 
exponent depending on two variables 

ln(p 2 + 2pcosh[A] + 1) a — > aab 

Z3f = 2hr3 f ° r 03f: b - bab ^ 

and substitution rules that lead to locally isomorphic (or patch equivalent) 
chains, and 

In 63/ , . a — > aab . . 

Z3/ = ^3" fOT Q3f: b - bba ^ 



where ©3/ = /OCOsh(A) + ^/[2 sinh(r/2) sinh(A)] 2 + 1. For r = 0, these ex- 
pressions reduce to the corresponding ones of the Ising quantum chain f20| . 
Note also that z$f = z-$f for pure isotropic aperiodicity. 



The substitution rule 



1 a — > abb , . 

b -> ababbb (?3) 



14 



belongs to case 3 of the classification above. A set of four strings, 



{(bb), (abab) , (abba) , (babbba) } , (74) 
is sufficient to define a string substitution with substitution matrix 



M„ = 



(2 2 1 2\ 
2 10 1 
111 

\0 1 2 3/ 



(75) 



We have A_ = 0, hence aperiodic anisotropy is irrelevant, but since X s = 2 
and ui2 = In 2/ In 5 > isotropic aperiodicity will be relevant. 



6 Extensions to correlated random disorder 

Aperiodic order, as generated by substitution rules, represents a natural, but non- 
trivial extension of crystalline and quasicrystalline order. Structures with this type 
of long-range order are certainly physically reasonable (perhaps in contrast to hi- 
erarchical systems) but of course also show a number of quite special properties in 
comparison to random systems. These include rescaling symmetries and the strong 
repetitivity of local patches due to their self-similar structure. Moreover, the order- 
ing is deterministic by construction and leads to zero entropy density. However, 



as has been argued in |23|, |24| , the thermodynamical properties of (quantum) spin 
models should be unaffected by most of these special properties, but depend only on 
the nature of the fluctuations present in the system. In fact, also our RG formalism 
may be applied to a more general class of models and in particular does not depend 
on the exact self-similarity of the substitution chains. For simplicity, we concentrate 
on aperiodic anisotropy (respectively the Ising case) . Consider a chain of couplings 
chosen according to the following random substitution rule 

g: at — ► perm(wj) (76) 

where perm(u>i) denotes a random permutation of the letters in Wi. The class of 
chains generated this way is clearly neither deterministic nor strongly repetitive, in 
fact it is, almost surely not repetitive at all. What is more, its entropy density is 
positive. Indeed, the nature of the fermionic spectrum is completely changed by 
the introduction of randomness: whereas it is typically purely singular continuous 
with a characteristic gap-structure for substitution chains, all these gaps vanish in 
(numerical) spectra of random-substitutions. As the only property that remains un- 
changed, the total fluctuation of the mean (reduced) coupling still decomposes into 
a superposition of a finite number of fluctuation modes, implicitly given through 
the eigenvalues of the substitution matrix. This property is also the essential in- 
gredient for our RG procedure, which depends on the substitution matrix rather 
than on the detailed form of the substitution itself. Consequently, neither the RG 
flow nor the fixed point structure are affected by introducing randomness into the 
substitution rules. Note however that the scaling exponent z in the marginal case 
depends on g in more detail and we only obtain analytical upper and lower bounds 
(1 < z\ < z < Z2 < oo) for random substitutions here p2| . Numerical results show 
that z may indeed vary within this interval and does not converge to a well defined 
limit. 

In many respects, random substitution chains with relevant fluctuation modes 
(in particular those with wandering exponent u>2 = 1/2) behave very similar to 
uncorrelated random chains. However, a characteristic difference is that for the 
latter only the asymptotic growth of the mean fluctuations is controlled by the mean 
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deviation exponent u = 1/2, while fluctuations of any order (up to u — 1) may be 
present with a non-vanishing probability on every lengthscale. The most significant 
consequence is the off-critical Griffiths phase observed in random quantum chains 
(|, ||, but not in aperiodic models |ll|. Also for random substitution chains with 
exponentially many realizations, no Griffiths phase should be present since for any 
non-critical values of the coupling constants there is a finite maximal size for "locally 
critical" patches. 



7 Discussion 

We extended an exact real space renormalization approach, originally formulated 
for Ising quantum chains to aperiodic XY quantum chains. This way, relevance cri- 
teria for aperiodic modulations have been obtained analytically for a second class of 
models. As predicted by the Harris-Luck relevance criterion, the geometrical fluc- 
tuation exponent plays the key role in the determination of the critical behaviour. 
However, the fluctuation exponent 102 of the sequence of ratios of consecuting e x 
and e y couplings, which matters for the XY models, may differ from the wandering 
exponent aj g of the sequence of interactions directly, which had been the crucial 
quantity in the Ising case. As a consequence, the relevance of the same aperiodic 
ordering of the coupling constants in general will be different for Ising quantum 
chains and isotropic XX models. In particular, quasi-periodic disorder, generated 
by substitution rules compatible with the cut-and-project formalism, is irrelevant 
for Ising quantum chains and most other Ising spin systems, but marginal for XX 
or XY chains. On the other hand, aperiodic XY anisotropy leads to identical fluc- 
tuations and the same critical behaviour as in the Ising case. 

The analysis of the RG fixed point structure and renormalization flows in par- 
ticular indicates that there is no discrimination between weak and strong aperiodic 
disorder in these models. The validity of the perturbative Harris-Luck criterion is 
thus extended to the case of strong modulations. Open questions remain mainly 
for relevant aperiodic disorder. Here the RG flows to the strong coupling limit and 
the critical scaling behaviour of several ensemble averaged quantities is dominated 
by rare events. A comparison of the resulting "aperiodic ground states" to the so- 
called "random singlet phase" postulated for uncorrelated random chains |j| would 
be of interest. For the Ising quantum chains, a first step into that direction has 
been done in We have shown that the RG approach may be also applied to 

random substitutions and does not rely too much on special properties of deter- 
ministic aperiodic systems. Let us remark that - especially in the Ising case - the 
structure of the RG is rather simple and an extension to ensembles of uncorrelated 
random chains should be possible somehow. The crucial question is whether the 
atypical means of quantities like the critical correlation function will be accessible 
within a RG scheme which works in the fermion representation. 

The renormalization approach leads to an exact determination of the scaling 
exponent z of the mass gap for arbitrary two-letter substitution rules and quite 
general classes of n-letter substitutions. The critical exponents connected with the 
scaling of the spectrum at A — may be calculated exactly, like the zero temperature 
specific heat, the susceptibility in a vanishing magnetic field in z-direction, or the 
localization length of the aperiodic hopping model at half filling. We have given a 
number of quantitative results as examples, mainly for marginal aperiodicity. 

The exact results obtained for the aperiodic XY models should be of use for 
the analysis (analytical and numerical) of more complicated aperiodic models, like 
Heisenberg and XYZ spin chains. As for analytical results, a natural step would be 
to introduce a non-vanishing transverse magnetic field in z direction. However, al- 
though this still leads to a free fermion model, a treatment by exact renormalization 
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as shown here at zero field does not seem to be a simple problem. 
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Appendix A: Spectra of string substitutions 

We are interested in the action of the transpose string substitution matrix M g on 
the vector (j, with 

/j, Si = A a # (si) + A 6 # 6 (s 4 ) + r(# ab - # 6o )0i) ■ (A.l) 

For fi = M*/x we find 

fi Si =A a # a (w Si ) + A h # b (w Si ) + r(# ab -# ba )(w Si ) . (A.2) 

Using 

#aK) = #aK)#«(*i) + #.W#l(«i) (A.3) 

#b{w Si ) = # b (w a )# a ( Si ) + # b (w b )# b { Si ) (A.4) 

we obtain the transformation rule ( |36| ) of the anisotropy parameters. 
On the other hand, since |u> a | is odd and \w b \ even, we may write 

- #ba)(w Si ) = (#ab ~ #ba)(w b ) ■ (2# 60 - (A.5) 

where #bo(si) gives the number of 6's in Sj = Subs^ with # a ( s n) even. Thus 
isotropic aperiodicity is clearly irrelevant if (# a & — #fc a )(wb) = 0. Otherwise, con- 
sider now the action of iVf' on fi. For A = we obtain 

[M*A] Sj = r(#„ b - #6.)K) • (2# 60 - #b)KJ (A.6) 

= K#a6 - #6a)K) ■ (2# 60 - ■ (2# fc0 - #*)(*) (A.7) 

and recognize A s = (2#b — #fc)( -!W h) as the desired RG eigenvalue. Since \w b \ is 
even, so is X s , isotropic aperiodic modulation is thus either irrelevant or relevant, 
but never marginal in this case. 
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